Perturbation Method for Classical Spinning Particle Motion: I. Kerr Space-Time 
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This paper presents an analytic perturbation approach to the dynamics of a classical spinning 
particle, according to the Mathisson-Papapetrou-Dixon (MPD) equations of motion, with a direct 
application to circular motion around a Kerr black hole. The formalism is established in terms of a 
power series expansion with respect to the particle's spin magnitude, where the particle's kinematic 
and dynamical degrees are expressed in a completely general form that can be constructed to infinite 
order in the expansion parameter. It is further shown that the particle's squared mass and spin 
magnitude can shift due to a classical analogue of radiative corrections that arise from spin-curvature 
coupling. Explicit expressions are determined for the case of circular motion near the event horizon a 
Kerr black hole, where the mass and spin shift contributions are dependent on the initial conditions 
of the particle's spin orientation. A preliminary analysis of the stability properties of the orbital 
motion in the Kerr background due to spin-curvature interactions is explored and briefly discussed. 

PACS numbers: 04.20.Cv, 04.25.-g, 04.70.Bw 



I. INTRODUCTION 

One of the earliest and on-going research interests in 
general relativity concerns the dynamics of extended bod- 
ies in the presence of strong gravitational backgrounds. 
Considering that virtually all astrophysical objects in the 
Universe, such as black holes, neutron stars, and other 
isolated massive bodies, have at least some spin angu- 
lar momentum in their formation, it is not difficult to 
surmise that an in-depth study of moving relativistic sys- 
tems with spin is a useful endeavour. A relevant example 
concerns the motion of rapidly rotating neutron stars in 
circular orbit around supermassive black holes like ones 
believed to exist in the centre of galaxies, which serve 
as candidate sources for emitting low-frequency gravita- 
tional wave radiation that may be detected by the space- 
based LISA gravitational wave observatory [1] . 

A first attempt to understand the dynamics of ex- 
tended bodies in curved space-time was put forward by 
Mathisson [2], who showed the existence of an interac- 
tion term involving the direct coupling of particle spin 
to the Riemann curvature tensor generated by a back- 
ground source. Steady progress was made since this 
first attempt, with a notable contribution made several 
years afterwards by Papapetrou [3], who proposed that 
the spinning particle exists within a space-time world 
tube containing its centre-of-mass worldline, where its 
associated matter field has compact support. In addi- 
tion, multipole moment contributions, i.e. beyond the 
mass monopole and spin dipole, to the extended objects 
full equations of motion were considered by Tulczyjew 
[4] and others, ultimately leading to the expressions ob- 
tained by Dixon [5, 6], with a self-consistent description 
for all multipole moment contributions to infinite order. 
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While the various theories of extended body motion in 
curved space-time differ with respect to the higher-order 
multipole moments, all of them recover the "pole-dipole 
approximation" identified initially by Mathisson and Pa- 
papetrou, which are satisfactory for most practical cal- 
culations, so long as the dimensions of the spinning body 
are small when compared to the background space-time's 
local radius of curvature. These truncated expressions of 
the full equations of motion are commonly known as the 
Mathisson-Papapetrou-Dixon (MPD) equations. 

There has been widespread interest in applying the 
MPD equations to the dynamics of classical spinning par- 
ticles in orbit around rotating black holes, as described 
by the Kerr metric [7-11]. In many ways, the Kerr back- 
ground is an ideal testing ground for the MPD equations, 
since both mass sources are spinning, which introduce in- 
teresting spin-curvature effects that impact upon the or- 
biting particle's overall evolution. Furthermore, it lends 
itself well to numerical simulations of deterministic chaos 
under extreme conditions [11-14], as well as studies of 
gravitational wave generation [15, 16] arising from spin- 
induced deviations away from geodesic motion. 

More formal study of the MPD equations have also oc- 
curred in various forms [17-19] , including a recent pertur- 
bative approach developed by Chicone, Mashhoon, and 
Punsly (CMP) [20], with application to the study of ro- 
tating plasma clumps propagating in astrophysical jets 
directed along a Kerr black hole's axis of symmetry. An- 
other application of the CMP approximation by Mash- 
hoon and Singh [21] determined analytic expressions for 
leading-order spin-curvature perturbations of a spinning 
particle's circular orbit around a Kerr black hole. This 
analysis is successful in reproducing the spinning parti- 
cle's kinematic behaviour compared to numerical simu- 
lations of the full MPD equations for situations where, 
for spin magnitude s and mass m, the M0ller radius 
[21, 22] for the spinning particle is s/m < 10~ 3 r, and 
r is the particle's radial distance away from the back- 
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ground mass source. However, this approximation starts 
to break down when s/ (mr) ~ 10~ 2 — 10 _1 for r = 10 M, 
where M is the Kerr black hole mass, suggesting that 
higher-order spin-curvature coupling terms are required 
to more completely describe the orbital motion. 



It was for this initial purpose that a generalization of 
the CMP approximation was very recently introduced 
by Singh [23] to incorporate higher-order analytic con- 
tributions to the perturbation approach for the MPD 
equations. This generalization has several nice features. 
For example, as a power series expansion with respect 
to the particle's spin magnitude, it can be extended to 
formally infinite order in the expansion. In addition, it 
leads to expressions that are background independent, 
and is fully applicable to arbitrary motion of the particle, 
without recourse to any space-time symmetries within 
the metric. As a result, this generalization is very ro- 
bust, with applicability for many distinct scenarios in 
theoretical astrophysics, such as the modelling of glob- 
ular clusters and other many-body dynamical systems 
in curved space-time, and also spinning particle interac- 
tions with gravitational waves, the results of which can 
be compared with existing treatments [24-26]. Further- 
more, this approach identifies the existence of a classical 
analogue for "radiative corrections" that shift the par- 
ticle's overall squared mass and spin magnitude due to 
higher-order spin-curvature contributions, a feature not 
thought about before. It would, therefore, be very useful 
to investigate the computational capacity of this general- 
ization when applied to circular motion in the Kerr back- 
ground. This is especially so in extreme conditions where 
a transition from stable to chaotic motion may be analyt- 
ically identified, for comparison with existing approaches 
[11-14] which use primarily numerical methods. 



The purpose of this paper is to present the generalized 
form of the CMP approximation for the MPD equations 
within the context of circular motion around a Kerr black 
hole, and explore the derived physical consequences. It 
begins with Sec. II, which displays the full MPD equa- 
tions, followed by a presentation of the formalism be- 
hind the generalized CMP approximation in Sec. III. Af- 
terwards, Sec. IV presents the formal application of the 
generalized CMP approximation to the case of circular 
motion around a Kerr black hole, up to second-order in 
the perturbation expansion parameter. This is followed, 
in Sec. V, by analysis of the predicted kinematic and 
dynamical properties of the perturbed system, includ- 
ing the predicted effective squared mass and spin magni- 
tude of the spinning particle. A general discussion of the 
main results obtained in this paper is found in Sec. VI, 
with a brief conclusion thereafter. The metric convention 
adopted is +2 signature with Riemann and Ricci tensor 
definitions following MTW [27], and geometric units of 
G = c = 1 are assumed throughout. 



II. MATHISSON-PAPAPETROU-DIXON (MPD) 
EQUATIONS 

The MPD equations of motion for the spinning par- 
ticle's linear four-momentum P^(r) and spin tensor 
S af \r) consist of 



DP* 1 
dr 

D g af 3 

dr 



va.fi 



l" S al3 . 



pa U P _ p/3 u a _ 



(la) 



(lb) 



where (la) describes the force applied due to spin- 
curvature coupling via S al3 (r) , the particle's four- velocity 
vector (r) = dcc M (r) / dr with affine parametrization t, 
and the Riemann curvature tensor R^ va p, while (lb) de- 
scribes the corresponding torque generated. As a result 
of (lb), the particle's four-momentum precesses around 
the centre-of-mass worldline. While it is possible to iden- 
tify r with proper time such that u^u^ = —1, it is use- 
ful to leave it unspecified at present. The differences 
between competing descriptions of extended objects in 
curved space-time arise from differing higher-order mul- 
tipole moment terms beyond the mass monopole and spin 
dipole moment, leading to additive contributions of the 
form and T a/3 [17, 21] in (la) and (lb), respectively. 
Specification of and T a/3 requires knowledge of the 
spinning object's energy-momentum tensor [5-7, 21], 
satisfying covariant conservation T^ v -. v = 0. For most 
practical purposes, however, the expression for (1) is suf- 
ficient. 

At present, the MPD equations as expressed in (1) 
are underdetermined, and require additional equations to 
completely specify the system. Following the approach 
from Dixon [5, 6], the orthogonality spin condition relat- 
ing the particle's linear and spin angular momenta ac- 
cording to 



S af3 P 8 = 



(2) 



is introduced, while the mass and spin parameters m and 
s naturally take the form 



m 



-P»P» 



2 <3jUI/ ^ 



(3a) 



(3b) 



It is very important to note that, while (3a) and (3b) are 
both technically functions of r, the MPD equations and 
(2) indicate that m and s are constants of the motion 
[20]. As well, the combination of (1) and (2) are known 
[9] to result in an expression for the four-velocity u M in 
terms of P^ and S af3 , such that 



u» = - 



P-u 



pn 



1 R via3 Pf S a(i 



2 m 2 + \ R aBpa Sp° 



, (4) 



where P ■ u is currently an undetermined scalar product 
relating the particle's internal clock with respect to r. It 
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becomes a self-evident confirmation from (4) that spin- 
curvature coupling displaces the particle's four-velocity 
away from a geodesic in curved space-time, leading to a 
dynamically rich interplay between the particle's centre- 
of-mass motion and its dynamical response due to spin- 
curvature interaction. The constraint equations (2)-(4) 
will prove very useful for ultimately deriving the gener- 
alized CMP approximation [23]. 



B. Generalization of the CMP Approximation 

1. Formalism 

The approach taken to generalize the CMP approxi- 
mation is to assume a power series expansion of the par- 
ticle's linear momentum and spin angular momentum, 
such that 



III. THE GENERALIZED CMP 
APPROXIMATION FOR THE MPD EQUATIONS 

A. CMP Approximation 



The approach taken by Chicone, Mashhoon, and Pun- 
sly in deriving the CMP approximation [20, 21] is to first 
assume that P M — mu^ = is a small quantity, where 
is the spin-curvature force. As well, the M0ller radius 
p [20-22] is chosen to be small, such that p — s/m <C r, 
where r is the distance from the particle to the back- 
ground gravitational source. This combination leads to 
the CMP approximation for the MPD equations, a series 
expansion to first-order in s, such that 



DP n 
dr 

DS af3 
dr 



0, 



(5a) 
(5b) 



where the spin tensor in (5b) is parallel transported 
within the approximation, and the spin condition (2) 
takes the form 



up 



(6) 



coinciding with the Pirani condition [28] relating the or- 
thogonality of the spin tensor to the four-velocity. 

The CMP approximation is a useful first step in an an- 
alytic perturbation approach to the MPD equations, with 
a remarkably accurate description for circular motion 
around a Kerr black hole [21] compared to the full MPD 
equations for s/(mr) <~ 10~ 3 and r = 10 M. However, 
it becomes clear that the CMP approximation breaks 
down as s/(mr) <~ 10~ 2 — 10 _1 for the same choice of 
r, which follows from the loss of torque information due 
to (5b), especially since the spin- induced modulation of 
the particle's r-dependent radial position found in the 
MPD equations is not present in the CMP approxima- 
tion. This weakness within (5) and (6) is suggestive of a 
more detailed and systematic approach that has resulted 
in the generalization to follow [23] . 



(3) 



(7a) 



3=0 

oo 



S»» (e) ee e £ ei S% = £ e? S^_ 1} , (7b) 

3=0 3=1 

where e is an expansion parameter to be associated with 
s, and and S^_^ arc the respective jth-order con- 
tributions of the linear momentum and spin angular mo- 
mentum in e. This implies that the zeroth-order expres- 
sions in e denote the dynamics of a spinless particle in 
geodesic motion. As well, the four-velocity is assumed to 
take the form 



S(e) ee ^ 



(8) 



3=0 



When substituting (7) and (8) into the MPD equations 
described by 

DP^(f\ 1 

— r-^ = --R" vafl u v {e)S a f ) {e), (9a) 

QT I 



DSa / {£) = 2 £ Pl«( £ ) </](£), 

dr 



(9b) 



where an extra factor of e is introduced in (9b) for con- 
sistency, it follows that the jth-order expressions of the 
MPD equations are 



DP^ i - 7-1 



dr 



nc«/3 3-1 

Ub (3-1) 



dr 



u fi] 

k) u (k) ■ 



k=0 



Given = m , where 



= _p(0) pA» 



(10b) 



(11) 



m — 1 \i 1 (0) ' 

it can be shown that the zeroth-order term in e is 
dr 

while the respective first-order terms following (10) are 

DP m 1 

= --R\ a eW( 0) S$, (13a) 



(12) 



DS 



a/3 
(0) 



dr 



= 0. 



(13b) 
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which is the CMP approximation. 



2. Supplementary Equations 



A complete specification of (10) requires determining 
u^.s as a function of the linear and spin angular momen- 
tum expansion components. In turn, this requires use of 
the supplementary equations (2)-(4), which have impor- 
tant consequences within the formalism. It is straight- 
forward to show that the spin condition (2) in terms of 
(7a) and (7b) is 



p(0) qliV 



E p i fe) So- 



ft) 



3 > 1 (14) 



k=l 



for the (j+l)th-order contribution in e, where the first- 
order perturbation in e is 



p(0) qflV 

M J (0) 



0. 



(15) 



Given that the mass and spin magnitude parameters to 
and s described by (3) are dependent on P M and S^ v , 
which arc represented by (7a) and (7b), respectively, it 
is possible to identify a classical analogue of a bare mass 
to defined by (11) and a bare spin s , according to 



are dimensionless jth-order corrections to to 2 , and s 2 ,, re- 
spectively. Given that to 2 and s 2 are already shown to be 
constant within the exact set of MPD equations, it must 
also be true that to 2 and s 2 are individually constant for 
each order of e. 

For the remaining supplementary equation (4), the 
four-velocity described by 



u"(e) 



P ■ u 

m 2 (e) 



1 S^(e)R viaP P^{e)S a P{e) 

2 m 2 (e)A(e) 



(21a) 



A( £ ) = 1 + 



1 



4ni2(:) R^S^(e)S^(e), (21b) 



can be determined as a series expansion in e, upon spec- 
ifying the yet undetermined scalar product P ■ u. With 
the particularly useful choice of 



P ■ u = —m(s), 



(22) 



_ q(°) 

2 °(0) 



(16) it is straightforward to determine that 



in analogy with the radiative corrections identified with 
the bare mass and spin parameters in quantum field the- 
ory. In this way, the MPD equations in perturbative form 
yield total mass and spin magnitudes due to the sum of 
"radiative corrections" to to and s , such that 



u„(e)ii"(e) - -1 



1 



4to 6 ( £ ) A 2 (e) 
l + 0( £ 4 ), 



(23) 



m 



'(e) = mUl + ^Te 



3 nm . 



3 = 1 



(17) 



where 



s 2 (e) = e 2 sl 1 + £ £ ^1 



i=i 



(18) 



R»{s) = S» v {s)R via0 P~<{s)S<*P{s). (24) 



where 



1 



TO, = — - 



y pu-k) P 



TO, 



2 f 
k=0 

j 



(fe) 



1 J 

,2 D (fc) ' 



k=0 



(19) 
(20) 



To at least third-order in e, the outcome (23) from (22) 
justifies the identification of t as proper time for parame- 
terization of the particle's centre-of-mass worldline. It is 
also straightforward, though tedious, to show explicitly 
from substituting (7) and (17) into (21) that the spinning 
particle's four-velocity in general form is 



■5 



u"(e) 



j=o 



e J u 



U) 



Mo) 
m 



+ s 



pM 
Ml) 



+ e^<! — 
m 



pA» 
^(2) 



2 TOl Mi) 



3 -4 

4 m i 



pp 

(o) 



a/3 



cA 1 ^ p p7 C" 

2 m 3 °(0) fl "1"F(0) D (0) 



m 



P^ 1 - _ m 2 
M3) o '"l M 2 ) 



+ 



2ml 



Rv^olQ 



EcM" p7 c a P ^ 2 cM^ p7 c 

Ml-n) 2^ >-fc) Mfc) ~ O 1 (°) (°) 1 



Ln=0 



+ 0(e 4 ), 



(25) 



which also satisfies (23) to third order in e. Following 
(25), it can be verified that [23] 



C. Solving for Linear Momentum and Spin 
Angular Momentum Expansion Components 



dr 



Dm) 
dr 



(26) 



up to third-order in e, with the use of (10). 



3. Perturbations of the M0ller Radius 



A useful consideration within the generalized CMP ap- 
proximation is the perturbation expression for the M0llcr 
radius p = s/m, since this quantity is relevant for iden- 
tifying the strength of the spin-curvature interaction for 
particles in the Kerr background while in circular orbit 
[21]. From previous work on chaotic dynamics [11—14], 
there is a strong suggestion that perturbations in the 
M0llcr radius may give some insight into determining the 
precise conditions for a transition away from stable mo- 
tion. A straightforward calculation shows that 



m(e) 



so 

2 



l + e 



m-i 



+ e 



2 

+ 0(e 3 )}, 



1-2 - 2\ 1 - 

(s 2 - m 2 ) - - s 



i m i 



ft -3 ml) 
(27) 



where the second- and third-order contributions in e 
due to the "radiative corrections" formally shift the per- 
turbed terms away from p = s /m . However, the pre- 
cise nature of the shift from p to p requires determining 
(27) in terms of a specific background. 



1. Local Fermi Co-ordinate Frame 

Having now presented the generalized CMP approxi- 
mation, the next step is to determine the linear momen- 
tum and spin angular momentum series expansion com- 
ponents for (7). This can be accomplished iteratively 
by solving the first-order perturbations with respect to 
zeroth-order quantities, and then computing the higher- 
order terms in a systematic fashion. This approach 
becomes particularly straightforward upon framing the 
problem in terms of the tetrad formalism and Fermi nor- 
mal co-ordinates [21] , the latter of which has the property 
that the corresponding metric in the neighbourhood of a 
freely falling worldline is locally flat. The leading-order 
metric deviations in Fermi normal co-ordinates are then 
proportional to the projected Riemann curvature tensor 
in the Fermi frame evaluated on the worldline. 

To begin, consider an orthonormal tetrad frame A m q, 
with the orthogonality condition 



r} & = 9„uX^&X V 
satisfying parallel transport 



dr 



= 



(28) 



(29) 



with respect to the general space-time co-ordinates A M , 
where a are indices for the Fermi co-ordinates X a defined 
in the neighbourhood of the spinning particle on a locally 
flat tangent space. Furthermore, the Riemann curvature 
tensor in the Fermi frame is described by 



R 



a/3y5 



R 



0- 



(30) 



= uj^ in the usual fashion and mak- 



By identifying A P q 
ing use of (28) and the first-order spin condition, it fol- 
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lows naturally that 

J* = \»&lf 0) = m A%, (31a) 

S(o) = ^"i ^(o) ' (31b) 

where P^ = m 5 a ^ and 5^ is a constant-valued spatial 
antisymmetric tensor whose components are determined 
from initial conditions. 

2. Leading Perturbation of Linear Momentum and Spin 
Angular Momentum 



is tempting to set s\ = in analogy with m 2 = 0, this 
is not well justified considering that s\ only needs to be 
covariantly constant following (26), and not necessarily 
zero. 

Fortunately, s\ can be determined separately by di- 
rectly solving (35), following a variation of an approach 
presented earlier [23] . Making use of the spin condition 
constraint equation with (14) for j = 1, there exist four 
equations 

A»S$-B V = 0, (39) 

where 



Determining the first-order perturbation in e for the 
linear momentum is very straightforward. Given (13a) 

and (29), it is shown that DPfa/dr = A^ a (dPfa/dr), 

leading to 



dP, 



(i) 



— _ _ F r?" 
d T ~ 2 °« (°) ' 



(32) 



which can be integrated immediately with the final result 
of 



; i ! - ~\>ft I ( F ^i%)dr. (33) 



It is interesting to note that, when (33) is contracted with 
P^ , the first-order mass shift contribution is identically 



m{ = 



(34) 



for a general space-time background, leading to simplified 
expressions for (25) and (27). 

The corresponding expression for the spin tensor, in 
contrast to the linear momentum, is somewhat more com- 
plicated to determine. It is important to first note that, 
from (10b) for j = 2, 



dr 



= 0. 



(35) 



In terms of the tetrad projection, 

— A & A D (l) 

= 2A^ Q A^5g+A^A^^ ); (36) 
where it follows from (14) for j = 1 that 



s?L = -—py's 



(!) qi] 



TO 



(0) 



(37) 



As for the components S 1 ,^-, in (36), they can be formally 
determined from using (20) for j = 1, such that 



Ail = pA» 

- (°) 



D = _pM o(0) 



(40a) 
(40b) 



With (39), it is possible to algebraically solve for the sff 
components in terms of the purely spatial components 



such that 



c(l) _ 
°01 — 


1 

A°~ 


A 2 


°12 


-A 3 


c(l) 
D 31 


+ B 1 


_ 


1 

A°" 


A 3 


c(l) 
J 23 


-A 1 


c(l) 
D 12 


+ B 2 


_ 

°03 — 


1 

A°~ 


A 1 


J 31 


-A 2 


c(l) 
°23 


+ B 3 



(41a) 
(41b) 
(41c) 



The remaining three components can then be determined 
by solving (35) in covariant form for the spatial compo- 
nents, leading to 



DS. 



(i) 



(i) 



(lT <h . 2«? 0) I%S$ = 0. (42) 

Upon substituting (41) into (42), there now exists a 
first-order inhomogeneous matrix differential equation to 
solve, with components 



-'12 

dr 



+ §a«Sg ) (r) 



dSilHr) , 1 

c 



dr 



5i 2 (r) ; 

$23 (t) • 



dS£>(r) , 1 



(lT +§7 l ^«(r) = *3i(r) 



(43a) 
(43b) 
(43c) 



where the antisymmetric spatial tensors a 4J , /3 y , and 7 y 
may each be r-dependent, depending on the choice of 
metric, and the Sij(r) are due to (40b). 



3. Higher-Order Perturbation Terms 



^(i) ~ 1 1 ^ 



(38) 



This, however, leads to a difficulty, in that (36) is still 
dependent on a yet undetermined parameter s\. While it 



Proceeding beyond the leading-order perturbations of 
P^ and S a @ is very straightforward. For the linear mo- 
mentum, the second-order expression is 
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1 

too 



F pa p/3 qkl 
/§fc? (1) (0) 



F pa C7/3 
" 07/3 



dr 



^ / (^^«^) + S <^> '^M-^mP^sgdr, (44) 



where 



<*1> = k[^ T)dT 



(45) 



is the time-averaged jth-order correction to the squared 
spin magnitude. To determine the second-order spin ten- 
sor, it is first shown from (10b) for j = 3 that 



nee" 

^(2) 

dr 



^fCTI S(L Rcfja/3 P? 0) S,£ . (46) 



TO 3 (0) ^(0) 



When expressed in terms of the orthonormal tetrad, it 
follows that (46) can be solved easily to obtain 

Sft = / S^ F R m sgdT. (47) 



IV. CIRCULAR MOTION NEAR THE EVENT 
HORIZON OF A KERR BLACK HOLE 

With the generalized CMP approximation of the MPD 
equations established, it can now be applied to the con- 
crete example of a spinning particle in circular motion 
around a Kerr black hole in the equatorial plane near its 
event horizon. It would be interesting to identify the gen- 
eral interplay between black hole spin and the dynamical 
response from the spinning particle. Since it is known 
that an orbiting particle whose spin co-rotates with the 
black hole spin leads to a repulsive force, while one of op- 
posite spin to the black hole leads to an attractive force 
[20], it is particularly useful to determine the black hole's 
spin dependence on the kinematics and dynamics of the 
particle's orbit due to contributions beyond the CMP ap- 
proximation considered earlier [21]. 

To begin, consider a Kerr background in standard 
Boyer-Lindquist co-ordinates X^ = (t,r,6,<f>), described 
in terms of a black hole mass M and specific spin angular 
momentum a — J/M. Then, for constant radius r, the 
solution to the geodesic equation for circular motion is 
[21] 



t 



— (l + aO K )T, 



N sin 9 



where 9 — tt/2 on the equatorial plane, 



fix = 



(48a) 
(48b) 

(49) 



is the Keplerian frequency of the orbit and 



N 



1 h 2a K 



(50) 



is a normalization constant. At r = 0, the boundary con- 
ditions are chosen such that t = <f> = 0. Given that A^q = 
dX^/dr is the particle's four- velocity vector, and that 
the orthonormal tetrad frame satisfies DX^a/dr = 0, it 
is possible to determine unit gyro axes A^j to describe 
the particle's local spatial frame. It is straightforward to 
show that [21, 29] 



A"n 



1 + aQ 



K 



A"t = - 



N 
L 

r A 



0,0,- 



N sin (9, 
sin (Ok t) , A cos (Ok r) , 0, 
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2M 



N 
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+ a Ok 



N 



lx / l-2aO K + ^ 



2M a 2 



(51a) 

(51b) 
(51c) 

(51d) 

(52a) 
(52b) 
(52c) 



for the circular orbit's energy E and orbital angular mo- 
mentum L. It is useful to describe a by the dimensionless 
parameter a = a/r, such that — 1/4 < en < 1 to incor- 
porate black hole spin (— M < a < M) that is both co- 
rotating and counter-rotating with respect to the orbital 
direction, where ro+ = M for a = 1 and ro_ = 4M for 
a = —1/4, each corresponding to the innermost (photon) 
radius allowed [30]. All dimensional quantities are then 
described with respect to M, such that 



9M 



■3(1- A^ 2 ) 



(53) 
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is the orbital radius near the event horizon and N > 
denotes the separation away from the innermost circular 
orbit. The parameters (52) are now expressible in terms 
of 

E = l-2(r 2 Q^) +a(r{l K ) = NE, (54a) 
L = (rO K ) [l-2 Q (rfi K ) + a 2 ] 

= N — , (54b) 
r 

I 



where A = r 2 [f — 2 (r 2 Q 2 ,) + a 2 ] = r 2 A 2 is the known 
function defined in the Kerr metric [30]. 

The exact expressions for F R^ 0& jj are listed in Ap- 
pendix A. For the special case of ./V — > and 9 = ir/2 
considered in this paper, the dominant nonzero contri- 
butions of the Ricmann curvature tensor F Rp iC)& in the 
Fermi frame are 



Fi? 6ioi « [2A 2 + r 2 Q 2 K -a(2rn K -a)] cos 2 (fi K r) = -^333, (55a) 

n 2 

F Rom ~ [ 2 ^ 2 + ?- 2 ^K (2 r fi K -a)] sin (fi K r) cos (fi K r) = - F i? 12 - 23 , (55b) 

3 ri 2 a 

Fi? 6ii3 ~ ( r - a) cos (fi K r) = - F R„m » (55c) 

^0202 ~ ^§ [l-2a(2rfi K -a)] = - F i? i3 i 3 > (55d) 
3fi 2 ^4 

Fi? 02l2 ~ ( r fi K - a) sin (fi K r) = - F i?63i3 . ( 55e ) 

Fi? 6303 « [24 2 + r 2 fi&-a(2rfi K -a)] sin 2 ^) = - F R im . (55f) 



A. First-Order Perturbations in e 

With all the ground work for the application of the gen- 
eralized CMP approximation now complete, it is possible 
to begin computing the linear and higher-order pertur- 
bations of the linear momentum and spin tensor for a 
spinning point particle in circular orbit. In this context, 
the expansion parameter e is associated with the unper- 
turbed M0ller radius po = so/m in unit of r, such that 
so/(m r) < 1. 

To begin, recall from (31a) that (51a) that the unper- 
turbed four-momentum components are 

P (°o) = ^(l + *m K ), (56a) 

P(o) = 0, (56b) 

P( ) = 0, (56c) 

PL = (56d) 



Obtaining the first-order perturbation (CMP approxi- 
mation) is very straightforward. Given (31b) and not- 
ing that the spinning particle is initially located on the 
x-axis of the Cartesian frame, the initial spin orienta- 
tion (9, <p) for S^ Q K coincides with the standard definition 
for the spherical co-ordinates (9, <fi) with respect to the 
Cartesian frame's z-axis. Therefore, it follows that the 
projected spin tensor components are [21] 

SV n = so sin 9 cos <f> , (57a) 

£31 = -s cos (9, (57b) 

^(0) = s o sin 9 sin <f> , (57c) 

leading to 
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so 
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= 0., 



( «o 



1 toq ff /" s 

2 rN A 
m E 



N 



m r 
so 
mo r 



cos 8 , 

sin(£lK t + — 0) — sin(OK t — 6 — <fi) 



cos(^k r + 9 — <j>) — cos(^k t — 9 — (J 
sin(£lK t + — <j)) — sin(riK t — 9 — <f>) 
cos 9 . 

I 



(58a) 

(58b) 
(58c) 
(58d) 

(58e) 

(58f) 



It is interesting to note the appearance of a beat structure 
in the sinusoidal functions of (58) , due to the initial spin 
orientation angles. As for the first-order perturbation for 



the linear momentum, this follows naturally from (33), 
resulting in 



3tooLo fJo_\ ( r fi K )( r fi K _ a ) \cos(n K T + 9)+cos(n K T-9)-2 cos8 
2 N a \to r J L 

l 1 ^- (— ) (rn K )(rn K -a)[sm(Q K T + 9)+sm(n K T-e) , 



2 iV 2 y too t 

3 Toq A ( s 



2 r AT 2 \ too r 

cos(0 + 0) — cos(0 — 



(rO K ) (rO K - a) 
(rfi K ) (r- - a) 



sin(fiK t + 9) + sin(fiK t — 
(0 K r + - <j>) - cos(0 K t - 6 - <j>) 



cos 



3 TOoSo (_£o\ / 0r x , ^ k _ a s L 0S m K T + m + cos (^ K r _ m _ 2 cos 
2 r N A \m r J L J 



(59a) 
(59b) 

(59c) 
(59d) 



which also exhibits a beat structure similar to what is 
found in (58). It is also interesting to note the relation- 
ship between the azimuthal and time component of 
in the form 



P o 

Mi) 



(r) 



E 

T 



(60) 



in agreement with the same computation performed ear- 
lier [21]. At this point, it is important to understand the 
conditions for the collapse of (59) when (rflx — a) =0. 
This condition can appear when 



M . 



(61) 



which is theoretically possible to reach when a — M [30] . 
However, it seems unlikely that such a possibility would 
arise in a realistic astrophysical context. 



B. Higher-Order Perturbations in e 



Evaluation of the second-order perturbation quantities 
is also straightforward, though rather involved. There 
are, however, relatively compact expressions for higher- 
order perturbation quantities in e that are required to 
evaluate the perturbed M0ller radius (27), namely the 
"radiative corrections" to the squared mass and spin 
magnitudes (19) and (20). The first computation of in- 
terest is s 2 , which is third-order in e according to (18). 
This is achieved by solving the first-order inhomogeneous 
matrix differential equation (43) for S^\t), followed by 
the algebraic equation (41) for the remaining components 
S^{t), which yields s\ via (20). Before proceeding, it 
is useful to first introduce a convenient notation for beat 
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functions taking the form 

Qf(ni Ok r , n 2 6 , n 3 cj>) = cos(niD,K t + n 2 9 - n 3 4>) 



Qt(ni K t , n 2 9 , n 3 (j>) = sin(niO K r + n 2 9 - n 3 <fi) 

± sin(mf2K t — n 2 9 — n 3 <f>) , 

(62b) 

which appear frequently in subsequent expressions 
throughout this paper. For the Kerr metric, it is shown 
that the nonzero a y , /3 U , and 7^ for (43) are 



while 



± cos(tiiOk t — n 2 9 — n 3 <p) , 

(62a) <Mt) 



Mr) 



3m rL / s \ 2 (r 2 0|)(rO K -«) 



4 TV 2 A \m r J (l + arfi K ) 
x Q s -(n K r,2M), (64a) 
J 23 (r) = 0, (64b) 



3 m rL ( s \ 2 2 



8 TV 3 V to r 



(r 2 K )(rO K -a) 



Q 



23 



0* 



Q+(2 K t , 2 , 2 0) — sin(2 K r - 2 

(64c 



(63a) 
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iV 

7a 2 " (l + arJ] K ) ' 
A 2 

-— (rO K ) (l + arO K ) , (63b) Then the solutions to (41) and (43) arc 



16 TO ° r ^m r; iV 3 (l+arU K ) 
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8 \TOor / N z 



{g+(2 K t , 2 , 2 $) - Q+(0 , 2 , 2 0) - 2 sin(2 K r - 2 0) + sin(2 0)] } 



-16 m ° r 



so 



Lo (r K ) (r K - a) 



m r J N 2 A (l + «rfi K ) 



{(2 Okt)Q s -(O k t, 20, $) - [Q c -(Okt, 26,$)- Q"(0, 20,0)]} , 
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16 V toq v 



N 3 



(r K ) (r K - a) 



x {(2 O K r)Q-(O K r,20,0)- [Q"(O k t , 2 9 , j>) - Q~(0 , 2 , 0)] } , 



--t m r 2 ( — ^_ ) — I (r K ) (r K - a) 
16 \m r ) N A 



x |Q+ (2 K t , 2 , 2 0) - Q+(0 , 2 , 2 4>) - 2 cos(2 K r - 2 0) - cos(2 0)J } 



(65a) 



(65b) 



(65c) 



(65d) 



(65e) 



(65f) 



When combined with (58), it follows that the first-order where 
shift in the squared spin magnitude is 
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= 16 



3 L (rn K )(rn K - a) 



(1 + arn K ) 
Q+(2Q K T,39,24>)-Q+(2n K T,§,2 



(2Q k t) Qt(2Q K T ,36 ,2 0) - Q+(2fi K T ,6 ,2 0) 
-Q+(ft K T,30,20) 



+ Q+(fiKT,fl,20)-Q+(n K T,3fl,O) + Q+(fi K T,e,O) 



- 2 



Q+(O,30,20) -Q+(O,0,20) +2 cos(3 0) - costf] } . 



(67) 



r 



The time-averaged expression for (66) over a cycle de- 
fined by the Keplerian frequency is 



1 



so 
m r 



(si) 



(68) 



where 



s 2 (r) dr 



2n 

3 L (rO K ) (ril K - a) 
2 



(1 + arVL K ) 
3 cos(2 0) - 1 



sin 2 9 cos 6 



(69) 



It is important to note that both (67) and (69) are 
well-behaved for the full range of 9 and 0. However, 
there exists the possibility for singularities to appear in 
the limit as (1 + ar fix) — * 0. This occurs when 



2/3 



a < 



(70) 



For the allowed radii permitted in the Kerr metric for 
photon orbits [30], it is clear that r > r c for a < 0, so the 
relevant expressions considered here can never become 
singular. 



It is more straightforward to solve for the remaining 
contributions to the perturbed M0ller radius, m 2 and s 2 , 
which are fourth-order perturbations in e according to 
(17) and (18). The first nonzero "radiative correction" 
to the squared mass magnitude (17) is 



m\(T) 



9 A 2 

8 A 4 



so 
mo r 



(r 2 n 2 K )(rn K -af 
r, 20,20) -Q+(O,20,2. 



{q+(2 k - 

cos(2 n K t - 2 0) - cos(2 0) | , (71) 



- 2 



according to (19), where its time-averaged expression fol- 
lowing the definition given in (69) is 



x (2 cos 2 6 - 1). (72) 
It is clear that the expressions (71) and (72) for m 2 . have 
no co-ordinate singularities due to 6 and 0, nor are there 
any physical singularities for reasonable choices for r. 

As for s\, it follows from evaluating (20) for j = 2 that 



3A 2 ( s 



«1(r) - ( — ] (r 2 K ) (rft K -a)(l + arQ K ) <2 C + (40 K r , 3 6 , 2 0) - Q+(4 K r , , 2 « 

8 A 4 \mo r j L 

-Q+(3f2KT,3^,2 0) + Q+(3n K r,^,20) + 3g+(3 Q K t , 3 , 0) + 5 Q+(3 r* K r,0,O) 

+ Q+(fi K r,3 0,-2 0) -3 g+(r» K r, 3 0,0)- Q+(f7 K r, 0,-2 0) -5 g+(fi K r, 0,0) 

- Q+(O,30,20)+Q+(O,0,2< 



(73) 



whose time-averaged expression is 

(4) - ~ (^) {r^D (m K -a)(l + arn K ) 
x sin 2 6 cos 6 (2 cos 2 - l) . (74) 



Again, (73) and (74) for s 2 , indicate well-behaved func- 
tions for all choices of 9 and 0, with no possibility of 
encountering singularities of any kind. It is interesting 
to note, however, that s\ is linear in so/(mo r), the same 
order as found in (66) for s 2 , which is somewhat un- 
expected. Furthermore, for A — > 0, it appears at first 
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glance that s\ will dominate over s\. Though it is difficult 
to identify the source for these unusual features, this may 
be indicative of another classical analogy to nonrenormal- 
izability in quantum field theory, where higher-order per- 
turbation terms in the generalized CMP approximation 
may possibly contribute to all orders of the expansion for 
certain quantities. This is a matter which may require 
further study in the future. 



V. NUMERICAL ANALYSIS 

At this point, it is useful to consider some numerical 
analysis of the main expressions for this paper, which 
are found in Appendix B of this paper. The purpose 
behind this procedure is to get a visual sense for how in- 
creasing the order of the perturbation expansion in the 
generalized CMP approximation influences the predicted 
physical behaviour of the spinning particle in the Kerr 
background. Since a purely numerical approach to the 
MPD equations does not allow for the clear identifica- 
tion of dominant contributions to the particle's orbital 
motion, this treatment provides an opportunity to glean 
some insight as to where a correspondence between the 
two approaches may occur. 

For all plots presented, r = 6M and 6 = <j> = ir/4. 
Particular attention is given to understanding the gen- 
eral stability of the spinning particle's motion while in 
circular orbit around a Kerr black hole, especially given 
the "radiative corrections" of the squared mass and spin 
magnitudes denoted by (19) and (20), respectively. 

For the purposes of this paper, p, = s /(m r) = 1CP 2 
and p = 10 _1 are considered throughout, where mo = 
10~ 2 M. It so happens that, for the given choices of r 
and mo, it follows that 

so = (l0 2 ml . (75) 

This implies that a realistic spin of s < m§ for solar 
mass black holes and neutron stars [13, 31] orbiting su- 
permassive black holes requires that 

M < 10- 2 f • (76) 

This upper bound given by (76) indicates that the choices 
of /i = 10~ 2 and p = 1CP 1 correspond to unrealistically 
large values [13, 14] for s , and consistent with values 
chosen in previous work [11, 12] suggesting chaotic be- 
haviour for the MPD equations. Therefore, any chaotic 
phenomena reported in this paper occurs under condi- 
tions not expected to be realized in a realistic astrophys- 
ical setting. 

To begin, consider the magnitude for the particle's co- 
ordinate speed 



v(s) = jgaV^V^e), (77) 



where 



and the u^(e) are given by (25). Since it must be true 
that < v < 1, it follows that any violation of this range 
of validity reflects a breakdown of the formalism's appli- 
cability. An exploration of (77) is presented in Figure 1, 
to first- and second-order in s, for the special cases of 
co-rotating (a = M) and counter-rotating (o = —M) 
extreme Kerr black holes. It is important to note that 
while r = 6M corresponds to stable orbital motion for a 
spinless particle when a = M, this choice for r only leads 
to marginally bounded orbits when a — —M [30]. 

For Figs. 1(a) and 1(b), the choice ofso/(mor) = 10~ 2 
shows that the spin-curvature force acting on the parti- 
cle's motion is almost exclusively due to the expression 
to first-order in e. In particular, the overall motion is 
stable throughout the range considered, with a variation 
on the order of 10~ 3 for Fig. 1(a) and 10~ 2 for Fig. 1(b), 
where the 0(e 2 ) expression only yields a 2 x 10~ 3 in- 
crease at the end of the plot compared to the first-order 
contribution alone. However, when s /(m r) = 10 _1 , 
the situation changes dramatically for both cases of a, as 
illustrated by Figs. 1(c) and 1(d). This is because the 
second-order expression in s introduces a rapid increase 
in the co-ordinate speed that approaches the v = 1 upper 
bound. It is particularly evident to see this in Fig. 1(d), 
which formally exceeds v(t) = 1 for r > 2000M. Such an 
outcome for Figs. 1(c) and 1(d) is consistent with prior 
numerical analysis on orbital stability [11, 13] when con- 
sidering large initial spin magnitudes sq. 

For immediate comparison, it is useful to now consider 
the M0ller radius p(r) = (s/m)(r) given by (27) for the 
same set of initial conditions. The purpose of this anal- 
ysis is to determine whether a correlation exists between 
the kinematic effects in v(t) with the anticipated dynam- 
ical contributions due to the spin-curvature interaction in 
p(r). Figure 2 is a plot of the M0ller radius in units of 
So/mo, up to 0(e 3 ), for a = M and a = —M. Compari- 
son with Figure 1 suggests that such a correlation exists 
between the two plots. The expression to third-order in 
e, which includes the mass shift contribution m\ as well 
as the second-order spin shift term s 2 ., has the effect of 
shifting the range of oscillation downwards within the 
plots as compared to the second-order expression alone, 
which oscillates with growing amplitude about p = 1. 

According to Fig. 2(a) for s /(m r) = 10~ 2 and 
a — M , the amplitude for M0ller radius grows very slowly 
when compared to Fig. 2(b) for s /(m r) = 10~ 2 and 
a = —M. Both the growth of the amplitude and the 
downward shift of the plots become more pronounced 
when examining Figs. 2(c) and 2(d) for a = M and 
a = —M, respectively, when so/( TO o r ) = 10 _1 . In par- 
ticular, the amplitude becomes many times larger than 
Po = sq /ma = 1 for both the co- rotating and counter- 
rotating black hole cases, which suggests that a certain 
minimum value for p(r) must occur before instability of 
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the orbital motion appears. Though the M0ller radius is 
not apparently a geometric quantity that must necessar- 
ily be positive-valued, it is interesting to note that the 
rapid increase in v(t) roughly coincides with the condi- 
tion that p(r) < for each of the plots in Figure 2. This 
may be a useful criterion for helping to determine the oc- 
currence of instabilities in the spinning particle's orbital 
motion. 

A further consideration involving the M0ller radius is 
to examine its time-averaged value (p) — (s /m) as a func- 
tion of the initial spin orientation angles and <j>. This 
leads to three-dimensional plots described by Figures 3 
and 4 for a = M and a — —M, respectively, which dis- 
play expressions to both second- and third-order in e for 
(p), and where So/(mor) = 10 _1 . It is evident that all 
the plots reflect an even function symmetry with respect 
to <f> = it. According to Figs. 3(a) and 4(a), the 0(e 2 ) 
expressions lead to a non-trivial peak and valley struc- 
ture in (p), while Figs. 3(b) and 4(b) for third-order in 
e effectively removes some of the structure for the region 
defined by < 9 < ir, leaving two peaks to dominate. 
While the shapes of the three-dimensional plots are ef- 
fectively unchanged when comparing between Figures 3 
and 4, going from a — M to a = —M leads to a ten-fold 
increase in magnitude, suggesting as expected that the 
counter-rotating black hole for r = 6M leads to a greater 
likelihood for encountering instabilities within the spin- 
ning particle's orbit. 

It is useful to briefly examine the linear momentum 
components P m (t), given (33) and (44). Figures 5-7 
display the radial, polar, and azimuthal components of 
the linear momentum, while Figure 8 displays the ra- 
tio P 3 (t)/P°(t), such as that described to first-order 
in e according to (60). For the radial component cor- 
responding to so/(rrior) = 10~ 2 , Fig. 5(a) for a = M 
shows that the expression to second-order in e introduces 
a slight contraction in the amplitude of -P 1 (r) before 
expanding outwards. This behaviour is also present in 
Fig. 5(b) for a = —M, though the outward growth is 
more pronounced, the beginning of which roughly cor- 
responds with the increase in v(t) in Fig. 1(b) start- 
ing at t = 3000M. Not surprisingly, the 0(e 2 ) ex- 
pression becomes dominant in Figs. 5(c) and 5(d) when 
So/(mor) = 1CP 1 , which also corresponds with the re- 
spective increases in v(t), as found in Figs. 1(c) and 1(d). 

The polar component corresponding to so /(trior) = 
1CP 2 is described by Fig. 6(a) for a = M and Fig. 6(b) 
for a = —M, which indicate a slightly net positive mag- 
nitude in P 2 (t) due to the expression to second-order in 
e. This outcome is somewhat surprising, since this sug- 
gests that the spinning particle will permanently leave 
the equatorial plane under these conditions. However, 
this may be more reflective of the choice for r, which 
hovers around the minimum value allowable for stable 
circular orbits. Again, this outcome is more pronounced 
for the case of so/( m o r ) = 10 _1 , as shown in Figs. 6(c) 
and 6(d). 

For the azimuthal component corresponding to 



So/{rriQr) = 10~ 2 , Figure 7 behaves similarly to that 
of Figure 5, particularly where it concerns Figs. 7(a) and 
7(b) for a = M and a = —M, respectively. That is, 
the 0(e 2 ) expression indicates a slight contraction in the 
amplitude of P 3 (r) prior to an outward expansion. Con- 
sistent with previous plots, Figs. 7(c) and 7(d) show a 
dominant growth of the amplitude of P 3 (r) for the ex- 
pression to second-order in e and So/(mor) = 10 _1 . 

Concerning the ratio P (t)/P (t), this is presented in 
Figure 8 for the case of so/(m r) = 10 -2 , where Fig. 8(a) 
refers to a = M and Fig. 8(b) corresponds to a = —M. 
As expected, the first-order contribution in e leads to 
a constant ratio in r, consistent with (60). When the 
second-order contribution in e is added, the ratio exhibits 
a slight contraction followed by an outward expansion, 
consistent with Figs. 7(a) and 7(b) for P 3 (t). Nonethe- 
less, it appears that the basic ratio remains constant for 
changing r. 

Finally, to explore the numerical properties of the spin 
tensor due to the generalized CMP approximation, con- 
sider the example of S 02 (t), as presented in Figure 9, 
which shows expressions up to third-order in e. It seems 
evident that the expression to second-order in e is dom- 
inant, as the it is clear that the 0(e 3 ) expression has 
no discernable impact on the amplitude. According to 
Fig. 9(a) for so/( TO o r ) = 10~ 2 and a = M, the ampli- 
tude remains constant around S 02 (t) — 0, while Fig. 9(b) 
for a — —M shows a gradual growth in the amplitude. 
Consistent with previous plots, this effect becomes more 
pronounced in Figs. 9(c) and 9(d) for s /(m r) = 10 _1 . 



VI. CONCLUSION 

This paper outlines the generalization of an analytic 
perturbation approach to the Mathisson-Papapetrou- 
Dixon equations for a spinning point particle, first in- 
troduced by Chicone, Mashhoon, and Punsly, with an 
application to circular motion around a Kerr black hole. 
The formalism shows the existence of "radiative correc- 
tions" to the particle's squared mass and spin magnitudes 
due to spin-curvature interactions, represented in power 
series expansion form. In performing the analysis, it is 
possible to semi-analytically identify the emergence of 
instabilities during the particle's orbital motion, which 
serves as a basis for a more precise treatment in the fu- 
ture. 

One of the underlying goals of the formalism presented 
in this paper is to determine the perturbed orbit of the 
spinning particle according to the generalized CMP ap- 
proximation, following the approach taken earlier [21]. 
However, to do this properly requires a modification of 
the equations of motion to incorporate dissipative effects 
due to gravitational radiation, which have not yet been 
taken into account. Such a modification would most cer- 
tainly require evaluation of the Teukolsky equations for 
determining the radiation effects corresponding to an adi- 
abatic inspiral for the spinning particle's orbit. This is 
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a non-trivial exercise with both conceptual and techni- 
cal challenges to still overcome. Once this is better un- 
derstood, a determination of the perturbed orbit due to 
spin-curvature interactions will be considered in a future 
publication. For now, a second paper on the general- 
ized CMP approximation in the Vaidya background is 
forthcoming [32] as a companion piece to accompany and 
compare with this paper. 
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Given that the nonzero Ricmann curvature tensor com- 
ponents in standard Boycr-Lindquist co-ordinates are 
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where the nonzero components of the Riemann curvature tensor 

F R~ . . a in the Fermi frame are listed as follows: 
E = r 2 + a 2 cos 2 9, (A. 14) ^ 

A = r 2 + a 2 -2Mr, (A.15) 
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FIG. 1: Co-ordinate speed v(r) of the spinning particle while in circular orbit around a Kerr black hole, for r = 6M and 
8 — (j> — tt/4. It is clear from Fig. 1(a) that v(t) to second-order contribution in e has no significant impact on altering the 
particle's orbital speed for so/{mor) — 10 -2 , though 1(b) indicates a slight increase over time. In contrast, Figs. 1(c) and 1(d) 
show that an instability occurs as so/(mor) = 10 _1 when considering the second-order contribution of e in v(t). 
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FIG. 2: M0ller radius p(r) = (s/m)(r) for r = 6M and 9 = (p = n/4, in units of so /mo- Fig. 2(a) shows that while the 
higher-order contributions in e lead to a slowly increasing amplitude in p, with a moderate increase found in Fig. 2(b). As 
so/(mo r) = 10 _1 , the amplitude increase becomes more pronounced for Figs. 2(c) and 2(d), where the second- and third-order 
contributions in e are noticeable. 




(a) a = M , 0{e 2 ) 




FIG. 3: Three-dimensional plot of 
the time-averaged M0ller radius (p) = 
(s/m) as a function of 6 and cj> for 
r = 6M and a = M. Fig. 3(a) shows 
a complicated peak and valley struc- 
ture to (p) that simplifies somewhat 
in Fig. 3(b), with two peaks present. 



(b) a = M, 0(e 3 ) 





FIG. 4: Time-averaged M0ller radius 
as a function of 8 and <j> for r — 6M 
and a = —M. It is clear that while 
the peak structure is essentially un- 
changed when compared to Figure 3, 
the magnitude increases ten-fold when 
going from a = M to a — —M. 



(b) a = -M , 0(e 3 ) 
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FIG. 5: Radial component -P 1 (t) of the linear momentum for r = 6M and 9 = <f> — tt/4. Fig. 5(a) shows a slightly growing 
amplitude due to the second-order contribution in e for so/(mor) = 10~ 2 , with a more moderate growth in Fig. 5(b). The 
amplitude grows much more rapidly for Figs. 5(c) and 5(d) as so/(mor) = 10 _1 . 



21 




0.0009 r 
0.0008 r 
0.0007 \ 
0.0006 - 
0.0005 - 
0.0004 j 
0.0003 j 
0.0002 - 
0.0001 - 
0.0000 j 
-0.0001 r 
-0.0002 ^ 



1 



iiiiiiiiiiiwi 



4000 5000 6000 7000 8000 9000 



T [M}_ 

(c) so/(mor) = 10 




"0 200 400 600 800 1000 1200 1400 1600 1800 2000 



(d) s /(m r) = 10" 1 , a = —A/ 



FIG. 6: Polar component P 2 (t) of the linear momentum for r — 6Af and — <j> = 7r/4. The second-order contribution in e 
introduces a slight non-zero value in the net magnitude for Figs. 6(a) and 6(b) with so/(mor) — 10~ 2 , whose average slope 
becomes more pronounced for Figs. 6(c) and 6(d) as so/(mor) = 10 _1 . 
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FIG. 7: Azimuthal component P 3 {r) of the linear momentum for r — 6M and 6 = = 7r/4. The second-order contribution 
in £ introduces a slight non-zero change in the amplitude for Fig. 7(a) with so/(mor) = 10 -2 , while Fig. 7(b) shows a more 
moderate growth in amplitude. This growth becomes strongly unbounded for Figs. 7(c) and 7(d) as so/(mor) — 10" 1 . 
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FIG. 8: Ratio of P 3 (r) to P°(r) for r = 6M and 6 = 4> = vr/4. When considering the expression to second-order in e, it is 
evident from Figs. 8(a) and 8(b) that the higher-order contribution becomes gradually unbounded compared to the constant 
ratio given by the first-order contribution in e. 
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FIG. 9: The S° 2 (t) component of the spin tensor for r — 6Af and 9 = <j> = 7r/4. In Fig. 9(a), there is virtually no contribution 
of the second-order and third-order contributions in e to the amplitude, while there appears a slight growth in amplitude 
for Fig. 9(b) due to these contributions. In contrast, Figs. 9(c) and 9(d) show a noticeable increase in the amplitude as 
so/(mor) = 1CP 1 due to the higher-order contributions in e. 



